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Abstract 



We discuss the relic density of the hghtest of the supersymmetric par- 
ticles in view of new cosmological data, which favour the concept of an 
accelerating Universe with a non- vanishing cosmological constant. Recent 
astrophysical observations provide us with very precise values of the rel- 
evant cosmological parameters. Certain of these parameters have direct 
implications on particle physics, e.g., the value of matter density, which 
in conjunction with electroweak precision data put severe constraints on 
the supersymmetry breaking scale. In the context of the Constrained 
Minimal Supersymmetric Standard Model (CMSSM) such limits read as: 
Ml/2 ^ 300 GeV-340 GeV, mo ^ 80 GeV- 130 GeV. Within the context 
of the CMSSM a way to avoid these constraints is either to go to the large 
tan/5 and /i > region, or make f/j, the next to lightest supersymmetric 
particle (LSP), be almost degenerate in mass with LSP. 
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I. INTRODUCTION 



During the last few years the knowledge of the cosmological parameters has started 
entering an era of high precision with far reaching consequences not only for cosmology, 
but for particle physics as well. The cosmic microwave background temperature is ac- 
curately known, Tq = 2.7277 ± 0.002 '^K, the Hubble parameter is determined with a 
relatively small error, Hq = 65 ± 5 Km/sec/Mpc, the baryonic mass density is precisely 
determined by big-bang nucleosynthesis, ^shf) = 0.019 ± 0.001, while the determina- 
tion of the age of the Universe from the oldest stars, as well as other sources, yield 
% = 14 ± 1.5 Gyr [1]. Very recent observations of type la supernovae (SNIa), as well 
as measurements of the anisotropy of the Cosmic Background Radiation (CBR) provide 
additional information favouring an almost flat and accelerating Universe, where the 
acceleration mainly is driven by a non-vanishing cosmological constant [2-4]. 

There is a growing consensus that the anisotropy of the CBR offers the best way to 
determine the curvature of the Universe and hence the total matter-energy density Qq 
[1]. The data are consistent with a fiat Universe, Qq = 1.0 ± 0.2, while we are confident 
that the radiation component of the matter-energy density, that is the contribution from 
CBR and/or ultra relativistic neutrinos, is very small [1,2]. Therefore the present matter- 
energy density can be decomposed principally into matter density flu and vacuum energy 

no = ^M + ^A- (1) 

There is also supporting evidence, coming from many independent astrophysical ob- 
servations, that the matter density weighs Um = 0.4 ±0.1 (see for instance Ref. [1] and 

references therein). The Qm, values are then restricted by the age of the Universe 
and by the value of the Hubble parameter through 

^""^Wol \/ n^(l - y) + nj,(y^ - y) + y ■ 

The constraints stemming from Eq. (2) are however less restrictive than those coming 
from the supernovae SNIa data. Recently two groups, the Supernova Cosmology Project 
[3] and the High-z Supernova Search Team [4], using different methods of analysis, each 
found evidence for accelerated expansion, driven by a vacuum energy contribution: 
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4 11 
^^A = 3^^^ + 3 ^ 6 ■ 

So, for Om = 0.4 ±0.1 this relation implies that the vacuum energy is non- vanishing, 
fl\ = 0.85 ± 0.2, a value which is compatible with a fiat Universe, as the anisotropy of 
CBR measurements indicate. Taking into account the fact that the baryonic contribution 
to the matter density is small, = 0.05 ± 0.005, the values for matter energy density 
flu result to a Cold Dark Matter (CDM) density Ocdm — 0.35 ± 0.1, which combined 
with more recent measurements [1,5] of the scaled Hubble parameter ho — 0.65 ± 0.05, 
result to small CDM rehc densities: 

OcDM V ^ 0.15 ±0.07. (4) 

From measurements of the ratio of the baryonic to total mass in rich clusters, smaller 
values for the mass density are obtained. This ratio is found to be Qb/^m ~ 0.15 [6,7] 
which entails to even tighter hmits ncouho^ = 0.12 ± 0.04 [8]. 

Such stringent bounds for the CDM relic density affect supersymmetric predictions 
and may lower the limits of the effective supersymmetry breaking scale, and hence the 
masses of the supersymmetric particles. In Ref. [9] within the framework of the string 
inspired no-scale SU (5) x C/(l) supergravity model, by relaxing the cosmological constant, 
regions of the parameter space compatible with Qcdm ^o^ « 0.2 were delineated, and 
phenomenological predictions for the sparticle spectrum were given. The relevance of the 
high precision cosmology to constrained supersymmetry was addressed in Refs. [10,11]. 
More recently the CDM relic abundance with non- vanishing cosmological constant, in the 
framework of the Minimal Supersymmetric Standard Model (MSSM), was shown to put 
limits on supersymmetric mass spectrum [12]. In fact it was shown that gauginos can be 
within LHC reach, if the recent cosmological data are used. As stated in Ref. [13] it is 
worth pointing out that while electroweak (EW) precision data are in perfect agreement 
with Standard Model (SM) predictions, and hence in agreement with supersymmetric 
models which are characterized by a large supersymmetry breaking scale Msusy [14], 
the data on ficDM ^o^ push Msusy to the opposite direction preferring small values of 
MsusY- Therefore EW precision data may be hard to reconcile with the assumption that 
the lightest supersymmetric particle (LSP or x), is a candidate for CDM [13]. 
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The method to calculate the rehc abundance of a Dark Matter (DM) candidate par- 
ticle in the Universe is outlined in Ref. [15]. In i?-parity conserving supcrsymmetric 
theories the LSP may be a neutralino, which is a good candidate to play the role of DM 
[16]. Many authors [17-30,10-12,31,32,35,36] have since calculated the rehc neutralino 
density. In the early works, only the most important neutralino annihilation channels 
were considered, but later works [27,28] included all annihilation channels. Also more 
refined calculations of thermal averages of cross sections were employed, which took into 
account threshold effects and integration over Breit-Wigner poles [33,34]. 

Our study in this paper is based on the CMSSM, which is motivated by Super- 
gravity, assuming universal boundary conditions for the soft supersymmetry breaking 
parameters, and in which the EW symmetry is radiatively broken [38]. Our strategy 
of calculating the neutralino relic density follows three steps: First the SUSY particle 
spectrum and the relevant couplings arc generated, according to the supcrsymmetric sce- 
nario mentioned above. Then the thermally averaged cross sections {av) are calculated 
in their non-relativistic limit, using analytic expressions. Finally we numerically solve 
the Boltzmann equation, which governs the evolution of the neutralino relic density, by 
using very accurate routines able to handle stiff problems of differential equations. Re- 
garding the calculation of the relic density, we solve the Boltzmann equation numerically 
by finding a proper boundary condition along the lines described in Ref. [26]. This is 
reminiscent of the WKB approximation; it yields very accurate results and differs from 
the standard approaches used in most works. We want to emphasize that for the sake 
of the effectiveness of our computational code we have chosen to use analytic results in 
order to calculate the amplitudes of the processes contributing to thermally averaged 
cross section {av) [27]. The price one pays, is that these analytic results break down 
in the vicinity of the poles or thresholds of the cross section. However the comparison 
of our results with those of other studies [24,28], which treat the problem of poles and 
thresholds in a more accurate manner by calculating numerically the thermally averaged 
cross section [33,34], shows that they are in striking agreement. This occurs, at least, in 
regions of the parameter space of the CMSSM where this comparison is feasible. 

The effect of the coannihilation between the LSP and the next-to-lightest super- 
symmetric particle (NLSP) is quite important and should be duly taken into account 
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[37,33,22,30,10]. The importance of coannihilation of the hghtest of the neutrahnos x, 
which in most of the parameter space of the CMSSM is a bino, with tr has been pointed 
out in Refs. [31,39]. X ^ coannihilation are of relevance for values of the parameters 
near the edge where x '^R ^"^^ almost degenerate in mass. In such regions of the 
parameter space the results reached using the ordinary methods, in which these effects 
are neglected, have to be properly modified to correctly account for the effect of the 
coannihilation. 

As a preview of our results: 
We have found that within the context of the CMSSM the recent cosmological data, in 
combination with EW precision measurements, lead to rather tight limits for the relevant 
supersymmetric breaking parameters mo, A'/i/2, provided the next to the LSP particle 
(tr) is not nearly degenerate in mass with the LSP. In this regime the only option to 
avoid these limits is to move to the large tan/? region, where acceptable relic densities 
can be obtained if the pscudoscalar Higgs mass is approaching twice the mass of the 
LSP. This case is consistent with 6 — > 57 and may be of relevance for models in which 
Yukawa coupling unification is enforced. 

In regions of the parameter space in which trS mass is close to that of the LSP, 
where coannihilation processes need be taken into account for the calculation of the 
actual neutralino relic abundance, such limits can be evaded. 

This paper is organized as follows: 
In the first section we give the basic formalism and discuss various details of our cal- 
culations. In section II and III we discuss the methodology we follow in solving the 
Boltzmann equation and give details of our numerical computation. In section IV our 
results for the LSP relic density are presented and regions of the parameter space con- 
sistent with the new astrophysical data are delineated. Towards the end of this section 
a discussion is devoted to the coannihilation effects. Finally we end up with the con- 
clusions. To facilitate the reader the supersymmetric conventions used throughout this 
paper are presented in the Appendix. 
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II. SUPERSYMMETRIC RELIC DENSITY 



Our aim is to calculate the cosmological relic density of the hghtest of the supersym- 
metric particles, which will be denoted by x throughout this paper. This we assume is 
one of the four neutralinos states. In supcrsymmetric models with i?-parity conserva- 
tion this particle is stable. The cosmological constraints on ^cdm discussed previously 
may impose stringent constraints on its mass, as well as on the masses of other super- 
symmetric particles which are exchanged in graphs, contributing to pair annihilation 
reactions 

constraining the predictions of supersymmetry. 

The basic ingredient in calculating the LSP rehc abundance is the calculation of the 
thermally averaged cross sections (avrei) for the annihilation processes xx ^ XY, which 
enter into the Boltzmann transport equation whose solution yields the mass density of 
the X particles at present epoch^. Vrei denotes the relative velocity of the two annihilating 
X's. Although these issues have been covered in numerous articles we will briefly repeat 
them from this stand too, in order to pave the ground for the discussion in the remainder 
of this paper. 

Our principal objective is to calculate the present LSP mass density 

Px = ^(^o) , (5) 

where Tq 2.7 °K is today's Universe temperature. This determines the LSP energy 
density fl^ = where pcrit is the critical density of Universe, is calculated by 
solving the Boltzmann equation given by 



where x = T/rriy^ and 



g = X{x) (q^ - ql) , (6) 



^ ~ T^h{T) ' ^° ~ T^h{T) ■ ^ ' 



We neglect at this stage slepton-x coannihilations and slepton-slepton annihilations. 
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In the equation above n denotes the number density of x's and no their density in 
thermal equihbrium. The latter is given by 



ksp^n T3 U^ju^ - 1) 

whose low temperature expansion (low x — T/m^) is 

In the equations above kspm is the number of the spin degrees of freedom. The function 
h(T) counts the effective entropy degrees of freedom, determining the entropy density of 
the Universe 

27r2 



- 45 h{T) , (10) 

which along with the effective energy degrees of freedom g{T), which determines the 
energy density 

p=^r^(r), (11) 

enter into the prefactor A(x) appearing on the right hand side of Eq. (6) 

X{x) ^ (^^G^j (h{T) + "^h'iT)) {avrei) . (12) 

Depending on the temperature T the content of the particles in equilibrium is dif- 
ferent. In our analyses we use the expressions for g(T), h(T) as given in Ref. [26]. In 
the region 40 MeV < T < 2.5 GeV, where the quark-hadron phase transition takes 
place, the values used for g{T), h{T) are those corresponding to a critical temperature 
Tc = 150 MeV as given in Ref. [18]. For a critical temperature = 400 MeV, also quoted 
in Ref. [18], we did not observe a substantial change in our final results concerning the 
LSP relic density. Recent lattice QCD results indicate that a first order phase transition 
takes place during the hadronization [40] . Using the corresponding data for the energy 
and entropy densities [41], no significant change is observed in our final results, as it has 
been also noticed in Ref. [18]. 

We postpone for later the details of the numerical scheme employed to solving the 
Boltzmann equation (6) and pass to discuss the thermal averages {aVrei) for the various 
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processes involved. At this point we follow Ref. [27] and express the non-relativistic 
cross sections for the annihilation processes % % — > X y in terms of helicity amplitudes 
as follows 

4 Svrs bf 

xJ2 (\A\'So)\' + k\A\'Po)\' + \A\'Pi)\' + \A\'P2)\')) , (13) 

h ^ ' 

where v is the relative velocity Vrei- In Eq. (13) the amplitudes J^i^^^'^Lj) depends on 
the hehcities of the final products denoted collectively by "/i" and the total cross section 
is obtained as an incoherent sum over the final hehcity states. The cross section will 
be expanded up to 0{v^^ terms and for this reason only S and P waves in the initial 
state are of relevance. The statistical factor Sf appearing in the denominator in Eq. (13) 
equals to 2! when the final particles are identical. The kinematical factor /?/ is given by 

-^,= {l- + + 2(m^-m^) y/^ ^ ^^^^ 

where s is the center of mass (CM) energy squared. 

Although our analysis in many respects resembles that pursued in Ref. [26] it differs in 
the particular method employed to calculate the thermal averaged cross sections, where 
we follow closely Ref. [27]. The results of the two approaches ought to be identical if it 
were not for the fact that some interference terms between graphs in processes involving 
Higgs particles in the final state or one Higgs and a Z-boson, were omitted. In our 
approach these terms are implicitly included in Eq. (13). 

Since the r.h.s. of Eq. (13) will be expanded up to terms O(v^), we need cast the 
helicity amplitudes into the following forms: 

A'^(i5o) = a[; + ay + ... (15) 
^'('Po,i,2)-&'(Po,i> + ... (16) 

The ellipses in the equations above include higher in v terms. 

Besides this the kinematical factor /3j has to be expanded, and also the CM energy 
squared variable s should be expressed in terms of the relative velocity v as given below 

= /3o + + 0{v') , s-^ = (1 - ^)/4m| . (17) 
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By using these, the cross section of Eq. (13) can be brought into the form 

vcr^a + ^v'^ (18) 
6 

with the constants a, b defined by the following expressions 

a^kJ2 f^o \ao\^ (19) 

h 

b^6kJ2{ l«ol' (A - §) + /^o (ao'*at + h.c.) 

h 

+ y ( \b\Po)\' + \b\Pi)\' + \b\P2)\' ) ) . (20) 

The prefactor k appearing in the equations above is given by 

k~^ - 128 TT 5/ m| . 

It is well known that the expansion in the relative velocity v breaks down near 
thresholds or poles. Concerning the kinematical factor ^9/ we write 



(21) 



where 

V4 mx my {mx + myf , . 

6 = , e = 1 . (22) 

mx + my 4 m^'^ 

This expansion obviously breaks down when e gets small, or equivalently when we are 
near the threshold 

2m^ = nix + . (23) 

Also singular are the expansions (19) and (20) when we are near an s-channel pole of 
a particle of mass mi into which xx fused to. The intermediate particle's propagator 
in this case is expanded as 

' <■ • (24) 



s-mj + imiFi m^^ 4:- Rj + iGj \ 4:-Rj + iGiJ ' 
where 

R:^^,G:^^. (25) 
m^ m\ 
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The expansion (24) holds as long as we are away from poles, otherwise the coefficient 
of the relative velocity squared gets large. The largeness of this factor is dictated by the 
narrowness of the resonance and the heaviness of the LSP. For the Z-boson resonance 
for instance, the corresponding rescaled width Gz is Gz ~ (230/m|) GeV^, which for 

?a 100 GeV yields Gz ~ 2.3 x 10~^ invahdating the expansion (24) on the resonance. 

Therefore near poles 



■m<, 



Mi 



(26) 



as well as near threshold, more sophisticated methods should be used, as those found in 
Rcfs. [33,34], for the non-relativistic expansion of the cross section in Eq. (18). We shall 
come back to this point later when discussing the LSP relic density. 

To make contact with the findings of Ref. [18] we write the cross section as 

1 



va — 



w{s) , 



(27) 



where Ei, E2 are the energies of the initial particles and s the total CM energy squared. 
Eq. (27) leads, up to 0{x), to a thermal averaged cross section (for details see Ref. [18]) 
given by 



(va) = -\ 
mi 



Wo + -(-2 Wq + Wq) X 



(28) 



where x — T/m^ and 



Wo = w{so) , w'q 



Am~^ ( 



/ dw\ 



(so = 4 m-") . 



(29) 



so 



By comparing Eqs. (13) and (27) we can have 



a — 



Wo , b 



{w'o - Wo) , 



(30) 



which can be used to cast Eq. (28) into the form 



(wa) ~ a + (6 a) x . 

2 



(31) 
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III. SOLVING THE BOLTZMANN EQUATION 



The coefficients a and b, appearing in Eq. (31), are calculated for each process 

XX^XY (32) 

where x is the lightest supersymmetric particle which we assume is one of the four 
neutralinos as said in previous sections. At low temperatures the particles in the final 
state may include ordinary fermions, gauge bosons or Higgses. 



The freeze out temperature Tf usually occurs for values oi Xf = 




0.05 and 



hence we can solve the Boltzmann equation (6) in the regime x < Xq by knowing the 
value of q{x) at a properly chosen point Xq > Xf which is not much beyond Xf ^. For 
temperatures T corresponding to x < Xq contributions of sparticles other than the LSP 
to g{T), h{T) are negligible, relative to LSP, and can be safely ignored. The reason is 
that any sparticlc's mass rrii is larger than and hence the relative Boltzmann factors 
exp [~—^;--^] arc suppressed in the region x < xq ^ 0.1. Hence only the contribution 
of the LSP is kept in the effective energy and entropy degrees of freedom functions g{T) 
and h{T) respectively^. 

Also, as stated previously, in the annihilation process in Eq. (32) only non- 
supersymmetric particles are considered in the final state. Although this is obviously 
correct at zero relative velocity of the initial particles (at threshold), since x is the LSP, it 
may not be the case at finite temperatures when x are adequately thermalized acquiring 
kinetic energies sufficient to produce heavier sparticles. Therefore channels which are 
forbidden at zero relative velocity may be activated at temperatures T. In this work 
we will follow the standard treatment and ignore contributions of all channels which 
are forbidden at zero relative velocity. This is justified by the following argument. The 



^The choice of xq is related to the particular method employed for solving the Boltzmann 
equation to be discussed later in this chapter. The resulting values of xq turn out to be 
around 0.1. 

^Obviously in regions where the coannihilation effects are important this approximation does 
not hold and the contributions of sparticles with masses close to mass of the x should be added 
to g{T), h{T). 



11 



values of x relevant for our calculation are x < 0.1 and as a consequence the correspond- 
ing temperatures are much smaller than m^. Therefore the initial state particles x are 

not adequately thermalized to activate a forbidden reaction. We can appeal to a more 
quantitative argument by recalling that in the forbidden region the thermally averaged 
cross sections are proportional to 

e-'^-/^ , (33) 

sec Rcf. [33], where yU_ depends on the masses of the final products X, Y . When for 
instance these have equal masses, say m2, and m2 > rriy^ this is given by = 1 — 
Therefore in the region a; < 0.1 the exponent in Eq. (33) drops rapidly, unless rrij^ is close 
to m2. This is what is intuitively expected; at low temperatures {T <^ m^) the thermal 
energies of the x's in the initial state are not sufficient to activate reactions in which 
the final products have masses well above their production threshold. Only when their 
masses are very close to threshold even a small amount energy is adequate to furnish 
enough kinetic energy to the initial particles to activate the reaction. On these grounds 
we therefore ignore the contributions of forbidden channels. This approximation is not 
expected to invalidate significantly our results. 

With this in mind the channels which contribute are (see also Refs. [26-28]), 

qq, II, W+W-, ZZ, ZH, Zh, ZA, W^H"^, HH, hh, Hh, AA, HA, hA, H+H' . 

q, I denote quarks and leptons, H, h, A denote the heavy, light and pseudoscalar Higgses 
respectively, while are the charged Higgses. The helicity amplitudes for the above 
processes have been calculated in Rcf. [27] as we have already discussed. Adjusting the 
results of that reference to conform with with our notation"^ we can calculate (aVrei)- 
Our numerical procedure then goes as follows: 

(i) Given the experimental inputs for SM fermion and gauge boson masses as well as 
couplings and supersymmetry breaking parameters, we first run two-loop Renor- 
malization Group Equations (RGE's) in order to define physical masses and cou- 
plings of all particles involved having as reference scale the physical Z-boson mass 
Mz. 



^Our notation differs slightly from that used in Ref. [27] (see Appendix). 
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(ii) We then calculate the coefficients a and b encountered in Eq. (31) for each of the 
processes mentioned before. 

(iii) We solve the Boltzmann equation to define the relic density at today's Universe 
temperature T ~ 2.7 °K . 

Regarding point (i) wc take as inputs the soft SUSY breaking parameters namely 
squark, slepton, Higgs soft masses, trilinear scalar couplings, gaugino masses as well as 
the parameters tan/? and sign(//). /i is the Higgsino mixing parameter. We assume 
CMSSM with universal boundary conditions at the unification scale Mqut- 

Although in our analysis we have enforced unification on gauge couplings at Mqut, 
the extracted values for the relic density are insensitive to this assumption and can cover 
cases where one abandons the naive gauge coupling unification scenario. In those cases 
the unification scale Mgut is defined as the point where gi and g2 meet. At this scale 
93 = fl'i,2(l + Ae). Ae ^ signals deviation from gauge coupling unification condition, 
which may be attributed to the appearance of high energy thresholds. Values of Ae 
of the order of 1% produce 5% variation in as{Mz), which however are not felt by the 
relic density. The reason is that the latter depends implicitly on as{Mz), through its 
dependence on sparticle masses, and therefore such small variations of ag have negligible 
effect on the relic density. Therefore our analysis can accommodate cases where one 
allows for small departures from schemes where gauge couplings unify at a common 
scale. 

Running two-loop RGE's for all couplings and masses involved, in the usual manner, 
we determine the parameters at the Z-pole mass which are necessary to calculate masses 
and couplings entering into the helicity amplitudes. Throughout radiative breaking of the 
EW symmetry is assumed. The magnitude of the fi parameter, but not its sign, as well 
as the Higgs mixing soft parameter mg are both determined at Mz via the minimization 
conditions of the one-loop corrected effective potential. 

All couplings and running masses are calculated in the DR scheme. Whenever needed 
these can be converted to their corresponding MS values. In a mass independent renor- 
malization scheme, as the DR, no theta functions enter into the RGE's to implement 
the decoupling of heavy sparticles at thresholds (see for instance Bagger et. al. in 
Ref. [14]). Therefore corrections to physical masses, which are calculated as the poles of 
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propagators, receive contribution from both light and heavy degrees of freedom. 

The pole masses of the third generation fermions are taken equal to = 175 GeV, 

M^"^^ = 5 GeV and Mf'^ = 1.777 GeV. From the pole masses wc can have the 
values of running masses at the pole, and then run the appropriate RGE's to have the 
corresponding running masses at the reference scale Mz- The b and r masses should 
evolve, according to the SU{3)c x U{l)em group, since M^''^^ M^°'^ are below Mz- Note 
that in the case of b and i-quarks, the two-loop QCD corrections relating pole and 
running mass are duly taken into account. In this way one obtains the values of the 
running masses at Mz, and from these the corresponding Yukawa couplings at the same 
scale in the DR scheme as demanded. 

Regarding Higgs boson masses, one-loop radiative correction to their masses are 
assumed through out this paper. The effect of the renormalization group improvement 
and leading two-loop corrections although important for an accurate determination of 
the Higgs masses does not significantly affect the values of the relic density. Only the 
location of the Higgs s-channel poles and the thresholds, whenever Higgses appear in 
the final state, are little affected. 

Radiative corrections to the couplings of the LSP to Higgses are not taken into account 
in this work. These can be important when LSP is a high purity Higgsino state [30], 
since a pure Higgsino state has no couphng to Higgs bosons. However in the CMSSM 
with universal boundary conditions for the soft masses at the Unification scale a high 
purity Higgsino state is hardly realized in view of negative results from SUSY particle 
searches, and the aforementioned corrections are not of relevance. 

We also assume that the LSP is the lightest of the neutralinos. At the stage (i) of 
collecting inputs for the calculation of the coefficients a, b we do not impose all existing 
experimental bounds on sparticle masses, especially those imposed on gluino and squark 
masses, some of which are conditional and model dependent. The reason for doing this 
relies on that we want to study the behaviour of the relic density ^^^q in as much 
enlarged portion of the parameter space as possible. Obviously the parameter space 
will shrink even more if additional experimental constraints are taken into account. We 
postpone a discussion concerning the experimental bounds used in our analysis for the 
following chapter. 
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Having all parameters at the scale Mz we pass to stage (ii) and calculate the coeffi- 
cients a and h (see Eqs. (19,20)) through which {aVrei) are calculated. As discussed in 
the previous section we have assumed non-relativistic approximation and have expanded 
up to 0{v'^) in the relative velocity v. However such an expansion breaks down near 
a threshold, or near a pole as discussed in the previous section. In order to quantify 
the notion of nearness to either a threshold or to a pole we first consider the threshold 
case. As is obvious from Eq. (22) we are on the threshold when the parameter e van- 
ishes, in which case the expansion of Eq. (21) breaks down. From this equation it is 
seen that the value of e signalling departure from the validity of the expansion in powers 
of V, is the one for which the coefficients of v'^ in Eq. (21) is unity. This occurs for 
eo = 0.125 which yields zq = — 0-94. Looking for a more reliable criterion we 

invoke Ref. [33] where results relying on more accurate analyses are compared against 
the standard schemes which wc arc using in this paper. From the figures displayed in 
the aforementioned reference we find that Zq ~ 0.95 not very far from the value quoted 
above. Therefore throughout our analysis we shall employ the following "near threshold" 
criterion 

Mx + My , , 

0.95 < — - < 1 . 34 
2m^ - ^ ' 

A similar analysis holds for the poles too. As is obvious from Eq. (24) the expansion 
breaks down when Rj = ^ is close to 2, unless the rescaled width Gj (see Eq. (24)) 
turns out to be large. The possible poles encountered are the Z, H, h and A resonances 
which have small rescaled widths unless the LSP is very light with mass 0(10) GeV. In 
our analysis we employ the following "near pole" criterion^ 

I 4 - I < 0.8 . (35) 

For values of the parameters leading to cither Eq. (34) or Eq. (35) the expansion of the 
cross section in powers of the relative velocity is untrustworthy and results based on such 
an expansion are unreliable. In those cases other more accurate methods should be used 
(see Refs. [33,34]). 



Outside this region the traditional series expansion, we use in this paper, and exact results 
are almost identical. See for instance Ref. [26]. 
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In the final stage (iii) we solve the Boltzmann equation (6). Knowing (aVrei) from the 
procedure outlined previously, and by calculating the functions g(T), h(T), h'(T), we 
can have the prefactor X{x) appearing in Eq. (6). At high temperatures, or same large 
values of a; = T/m^, above the freeze-out temperature, the function q{x) approaches its 
equilibrium value qo{x) (see Eq. (7)). A convenient and accurate method for solving the 
Boltzmann equation is the WKB approximation employed in Ref . [26] . This relies on the 
observation that X{x) is a rather large number of the order of 10^ or so, or even larger 
in some cases. Due to the largeness of X{x) an approximate solution is 

Obviously this holds for values of x for which is smaller than unity. In our numerical 
procedure we find a point xq for which 

^W.O.l. (37) 

For larger values of x, this ratio becomes even smaller while for smaller values increases 
rapidly invalidating the approximation (36). This rapid change of the aforementioned 
ratio is mainly due to A typical sample is shown in Table I where for an LSP mass 
fa 100 GeV, for a top mass equal to 175 GeV and for masses of the Higgses h, H, A, 
equal to 100, 250, 270 and 300 GeV respectively, we list its values for x in the 
range 0.03-0.08. One observes that ^ increases by almost 10 orders of magnitude from 
X — 0.08 down io x — 0.03 while X{x) remains almost constant in this interval. With a 
typical value of X{x) 10^° this ratio turns out to be ^ for values of x around 0.06. 

Given the point Xq, defined in Eq. (37), we numerically solve the Boltzmann equation, 
in order to obtain solutions in the region x < Xq, having as boundary condition 

g(xo) = 1.1 go(a^o) ■ (38) 

The omitted 0{1/X^) terms in Eq. (36) yield corrections which are less than 5%. There- 
fore this scheme yields very accurate results. 

The numerical solution is found by use of special routines found in IMSL FORTRAN 
library, which are eligible to handle stiff differential equations such as the Boltzmann 
equation. Similar routines can be presumably found in other libraries too. However 
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it is important to stress that the choice of the right routine and accuracy is of great 
importance. Due to the fact that q{x) varies rapidly for x < Xq a high degree of accuracy 
is demanded which makes other routines being either extremely slow or unable to reach 
convergence. 

To implement the numerical solution of the Boltzmann equation we need as inputs 
the function X{x), defined by Eq. (12), whose values are known provided g, h, h', as 
well as (avrei), are calculated. The effective number and entropy degrees of freedom 
functions g and h respectively are calculated in the way described in the first section. 
The thermal integrals needed for their calculation are found by invoking fast and rehable 
integration routines found in the same Fortran library IMSL. Their correctness has been 
checked by comparing our findings against those of other packages. The LSP mass and 
the masses of Higgscs and the remaining SM particles arc needed in order to calculate 
the aforementioned functions. Therefore we first run to get the values of all parame- 
ters involved at the physical scale Mz, as well as all physical masses among these the 
radiatively corrected Higgs masses. These inputs are also used in order to calculate all 
relevant sparticle couplings to other species necessary to calculate the coefficients a, b 
and hence (aVrei), for each one of the processes involved. 

After this short description of our numerical procedure we pass to discuss how this 
machinery is implemented to infer physics conclusions for the LSP relic density. 

IV. THE LSP RELIC DENSITY 

Following the numerical procedure outlined in the previous section we are ready to 
embark on discussing the predictions for the LSP relic density. As discussed in section 
I we have in mind minimal supersymmetry with universal boundary conditions at the 
unification scale for the soft breaking parameters and radiatively induced EW symmetry 
breaking. Therefore the arbitrary parameters are mo, M1/2, Aq and tan/5. The value 
of fi is determined from the minimization conditions of the one-loop corrected effective 
potential. These also determine the Higgs mixing parameter m|. The sign of /i is 
undetermined in this procedure and in our analysis both signs of /i are considered. 
Therefore in this scheme the /i value as well as m| are not inputs. 

At the first stage for each point in the parameter space we collect outputs including 
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all parameters relevant for the calculation of the relic density, such as couplings and 
physical masses, in the way prescribed in the previous section, without imposing any 
experimental constraints. However we certainly exclude points that are theoretically 
forbidden, such as those leading to breaking of lepton and/or color number, or points 
for which Landau poles are developed and so on. We also exclude points for which the 
LSP is not a neutralino. In subsequent runs the above inputs are used to determine the 
X relic density solving the Boltzmann equation as outlined in the previous section. 

In our analysis we should exclude points of the parameter space for which violation 
of the experimental bounds on sparticle masses is encountered. We use the bounds of 
Ref. [42] 



Neutralinos : m^o > 33 GeV 
Charginos : m~+ > 95 GeV 

Al 



Sleptons : ruf^ > 71 GeV , (m^o < 20 GeV) 

"^Afl ' "^efl > 84 , 89 GeV , (m/ij^^ga > m^o + 10) 
rriiy^ > 43 GeV 

Higgses -mho > 81 GeV, (light scalar) 

niA > 81 GeV , mH± > 69 GeV . (39) 

At this stage we do not exclude yet points which violate the gluino g and squark q mass 
bounds 

m-g > 173 GeV, mq > 176 GeV. (40) 

Then for each point of the parameter space for which the above experimental constraints 
are obeyed we calculate the x relic density. 

Prom our outputs we have found that the chargino bound is the most stringent of 
all listed above, in the parameter space examined. The only exception is the light Higgs 
boson mass, which outstrips the chargino bound for very low tan (3 values. The gluino 
and squark mass bounds quoted before, if subsequently imposed, are found to be weaker 
than the chargino bound. Only a small region of the parameter space which is allowed by 
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the chargino mass constraint is excluded when one enforces the bound rrif^ > 176 GeV 
on the hghtest of the top squarks. 

Before embarking to discuss our physics results we should stress that in our scheme 
we have not committed ourselves to any particular approximation concerning the masses 
or couplings of sectors which are rather involved, such as neutralinos for instance, which 
are crucial for our analysis. Therefore we do not only consider regions of the parameter 
space in which the LSP is either purely a B (bino) or purely a Higgsino, but also regions 
where in general the LSP happens to be an admixture of the four available degrees of 
freedom^. Regarding the LSP's mass we note that for large values of it many channels 
are open but for small values {m^ < 40 GeV) only channels with fermions, except the 
top quark, in the final state are contributing. In these processes the exchanged particles 
can be either a Z-boson and a Higgs in the s-channcl, as well as a sfermion / in the 
t-channel. Higgs exchanges arc suppressed by their small couplings to light fermions, and 
sfermion exchanges are suppressed when their masses are large. Then the only surviving 
term, for large values of squark and slepton masses is the Z-boson exchange. However 
in the parameter region where the LSP is a high purity bino, this is not coupled to 
the Z-boson resulting to very small cross sections enhancing dramatically the LSP relic 
density. Therefore in considerations in which the LSP is a hght bino^, large squark or 
slepton masses are inevitably excluded since they lead to large rehc densities. If one 
relaxes this assumption and considers regions of the parameter space in which the LSP 
is light but is not purely bino, heavy squarks or sleptons are in principle allowed. When 
LSP is light the only open channels are those involving light fermions in the final state. 
Then the annihilation of LSP's into neutrinos for instance, a channel which is always 
open, proceeds via the exchange of a Z-boson which is non-vanishing and dominates the 



^The case of a Higgsino-like LSP has been pursued in Refs. [30,32] where the dominant ra- 
diative corrections to neutralino masses are considered. Analogous corrections to couplings 
of Higgsino-like neutralinos to Z and Higgs bosons are important and can change the relic 
density by a factor 5 in regions of parameter space where LSP is a high purity Higgsino state 
[30]. However this case is not realized within the CMSSM with universal boundary conditions 
for the soft masses. 

^This happens when S> Mw, with Mi small ^ Mw- 
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reaction when mo is sufficiently large, due to the heaviness of sfermions. This puts a lower 
bound on (va{xx ^ f f)) and hence an upper bound on ^^hl which can be within 
the experimental limits quoted in the introduction. On these grounds one would expect 
that by increasing mo, while keeping M1/2 fixed and low, there are regions in which the 
relic density stays below its upper experimental limit. Although such corridors of low 
Ml/2 and large mo values^ are cosmologically acceptable they are ruled out by the recent 
bound put on the chargino mass. Hence the possibility of having a light LSP and a heavy 
sfermion spectrum is excluded. 

We have scanned the parameter space for values of mo, M1/2, Aq up to 1 TeV and 
tan/3 from around 1.8 to 40 for both positive and negative values of pL. The top quark 
mass is taken 175 GeV. In figure 1 we display representative outputs in the (M1/2, mo) 
plane for fixed values of Aq and tan f3. Both signs of the parameter fi arc considered. 
In the displayed figure Aq = and tan/? = 5, 20. The five different grey tone regions 
met as we move from bottom left to right up, correspond to regions in which /?.q takes 
values in the intervals 0.00 - 0.08, 0.08 - 0.22, 0.22 - 0.35, 0.35 - 0.60 and 0.60 - 1.00 
respectively^. In the blanc area covering the right up region, the relic density is found to 
be larger than unity. The boundary of the area excluded by chargino searches, designated 
by mg < 95 GeV, refers to the bound quoted in the beginning of this section. 

In these figures whenever a cross appears it designates that we are near either a pole 
or a threshold, according to the criteria given in the previous section. In these cases the 
non-relativistic expansions used are untrustworthy and no safe conclusions can be drawn. 
For low values of M1/2 crosses correspond to mainly poles, which are either a Z-boson or 
a light Higgs, while for higher values, where LSP is heavier and hence more channels are 
open, these correspond to thresholds. The grey area at the bottom labelled by "TH", 
which usually occurs for low values mo < 150 GeV, is excluded mainly because it includes 
points for which the LSP is not a neutralino. In a lesser extend some of these correspond 
to points which are theoretically excluded in the sense that either radiative breaking 



^These corridors of low M1/2 and large mo values have been also presented in Ref. [28]. 

^These regions are chosen in accord with new and old bounds on Qy. /iq which have been cited 
in the literature. 
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of the EW symmetry does not occur and/or other unwanted minima, breaking color or 
lepton number, are developed. Prom these figures it is seen that as tan /? increases from 
tan j3 = 5 to tan (3 = 20 the region for which the LSP is not a neutralino is enlarged. 
This is due to the fact that by increasing tan/? the stau sfermion tr becomes lighter, 
since its mass, as do the masses of all the third generations sfermions, depends rather 
strongly on tan (3 (and also on ^o)- Although not displayed, similar is the case when one 
increases the value of the parameter Aq. 

For fixed M1/2 > 150 GeV the rehc density increases, with increasing mo, due 

to the fact that cross sections involving sfermion exchanges decrease. Thus the area 
corresponding to ^x^o < 0-22 concentrates to the left bottom of the figure. In this 
region mo < 200 GeV. For fixed mo the rehc density il^ ^0 < 0-22 also decreases with 
increasing M1/2, since an increase in M1/2 enlarges squark and slepton masses as well 
yielding smaller cross sections. If M1/2 is further increased the LSP will eventually cease 
to be a neutralino. 

In figure 2 and for fixed values of the parameter Aq and tan /? we plot the LSP relic 
density as function of the soft scalar mass mo for values of M1/2 = 170, 200, 400 GeV 
respectively. The value M1/2 = 170 GeV has been chosen close to the lowest allowed 
by the recent chargino searches, and avoids poles or thresholds. It is obvious from this 
figure that for higher tan/3 values f^x^o S®^^ lowered, for fixed mo, leaving more room 
for larger mo and hence for sfermion masses. The abrupt stop in some of the displayed 
lines, towards their left endings, is due to the fact that the LSP ceases to be a neutralino 
for sufficiently low values of mo- 

In figures 3 and 4 we plot the LSP relic density as function of the parameters Aq and 
tan j3 respectively by keeping, in each case, the other parameters fixed. In figure 3 we see 
that for a relatively large value of the parameter mo = 200 GeV, and for all cases shown, 
the rehc density takes unacceptably large values. Although we have only depicted the 
case mo = 200 GeV this holds true even for larger values of mo, provided that M1/2 stays 
larger than about 150 GeV. 

The behaviour of f^x^O' ^ parameter tan/3 varies from 2 to 35, is depicted 
in figure 4. Keeping the parameters Aq, mo fixed one observes that for large values 
of tan/3 > 30, fl^h^ gets smaller falling below 0.22 even for large values of the soft 
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parameter Mi/2. The reason of getting small relic densities for such large values of tan (3 
is due to the fact that in these cases the pseudoscalar Higgs boson A has a mass close to 
2 m^, and thus its exchange dominates in the production of a fermion-antifermion pair 
in the final state. This, along with the fact that Arf and Abb vertices are proportional 
to tan/3, enhances the relevant cross sections, resulting to small relic densities within 
the allowed cosmological hmits. This behaviour agrees at least qualitatively with the 
findings of Ref . [27] (see figure 4 in that reference) . 

In figure 5 the LSP relic density is plotted as a function of the parameter M1/2 
for values of Aq, tan (3 shown on the figures, and for mo = 150 GeV (solid hne) and 
mo = 200 GeV (dashed line). The crosses denote points for which poles or thresholds 
are encountered. It is obvious in these figures the tendency for the LSP rehc density to 
increase as M1/2 increases especially for values M1/2 > 300 GcV. In this region, and 
for fixed Ml/2, we observe that flx^o decreases as tan/? is increased from tan/3 = 5 to 
tan/5 = 20. 

So far in our analysis we have not studied neutralino-stau coannihilation effects, 
which if included can lower the values of the neutralino relic density in some regions of 
the parameter space. However as we shall see, even in those cases our calculation of relic 
density can be used to estimate with fair accuracy the actual relic density by using the 
results of Ref. [31]. 

These coannihilation processes are of relevance for values of the parameters for which 
< rrif^ < 1.2m^, that is near the edge where x and ff> are almost degenerate in mass. 
Since so far in our analysis we have neglected such coannihilation effects, the conclusions 
reached are actually valid outside the stripe < nif^ < 1.2m^. Inside this band x — f 
coannihilations, and also / — / annihilations, dominate the cross sections, decreasing x 
relic densities leaving corridors of opportunity to high M1/2 and mo values as emphasized 
in other studies [31]. Thus depending on the inputs mo, M1/2, ^0, tan/3 and the sign of 
// we can distinguish two cases: 

(i) 1.25 rriy^ < mf^ and 

(ii) < rrifj^ < 1.25 m^, 

which are both compatible with having LSP as one of the neutralino states. In region 
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(ii) the stau tr is nearly degenerate in mass with x a-nd X ~ ^ coannihilation effects, 
and to a lesser extend f — f, e — e and fi — ji annihilations, play an important role [31]. 
We shall call this "coannihilation" region to be distincted from region (i) which will be 
designated hereafter as "coannihilation free" region. 

We shall first discuss the region (i) in which such effects are negligible and the ordinary 
way of calculating x i"elic densities, with the omission of the coannihilation processes, is 
very accurate and reliable. In the coannihilation free region upper limits on M1/2 and 
mo can be estabhshed by imposing the cosmological constraint 0.08 < JI^/iq < 0-22, 
which are more strict than those discussed so far. In fact within the coannihilation free 
region we find that for low and moderate tan j3 the upper bounds on these parameters are 
Ml 12 < 340 GeV, mo < 200 GeV. The upper limit set on mo is correlated to the value 
of Ml/2, and is almost insensitive to the value of the parameter Aq and tan/? as long as 
the latter does not get values larger than about ~ 30. For instance the upper bound 
~ 200 GeV on mo is reached when M1/2 ~ 140 GeV, the lowest allowed by chargino 
searches, but it is lowered to ~ 130 GeV when M1/2 — 340 GeV. This behaviour is 
very clearly seen in the scattered plots shown in figure 6. The sample consists of 4000 
random points that cover the most interesting part of the parameter space, which is 
within the limits: 1.8 < tan/3 < 40, 150 GeV < M1/2 < 1 TeV, |^o| < 500 GeV 
and mo < 500 GeV^°. Prom the given sample only points which lie entirely within the 
coannihilation free region are shown. Also points which lead to relics larger than 1.5 
are not displayed in the figure. The experimental bounds discussed before, restrict by 
about 40% the values of the allowed points. The points shown arc struck by a cross (x) 
when mo < 100 GeV, by a plus (+) when 100 GeV < mo < 200 GeV and by a diamond 
(o) when mo exceeds 200 GeV. It is obvious the tendency to have M1/2 < 340 GeV 
in the cosmologically interesting domain which lies in the stripe between the two lines 
at 0.08 and 0.22. Actually except for a few isolated cases, which correspond to large 
tan P as we shall see, all allowed points are accumulated to values M1/2 < 340 GeV and 



Higher values for mo are of relevance only for low Mx/2 values, already ruled out by the 
recent experimental bounds on chargino masses. Also since /iq does not depend strongly 
on Aq for Ml/2 > 1^0 GeV, as it can be realized from figure 3, it suffices to focus on values 
1^1 < 500 GeV. 
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mo < 200 GeV (crosses or pluses). 

As a side remark, we point out that the coannihilation free region under discussion 
overlaps with the color and charge breaking (CCB) free region as long as the parameter 
Ml/2 stays less than ~ 300 GeV [31]. 

Anticipating a forthcoming discussion on EW precision data, we designate the region 
of Ml/2, which is rather favoured by EW precision measurements. This is shaded in grey, 
which progressively becomes darker as we move to larger M1/2 values, where the SM 
limit is attained. In the coannihilation free region, these upper limits set on M1/2, mo 
can be only evaded when tan /3 takes large values ~ 30 and // is positive. Higher values 
of tan P can be also obtained at the expense of changing the input value for the bottom 
pole mass as we are discussing below. In the aforementioned cases the pseudoscalar 
Higgs A has a mass approaching 2 m^, and the Att, Abb couplings are large as being 
proportional to tan/?. Both effects make the pseudoscalar Higgs exchange to dominate 
the reactions X X ^ ^ ^ X X ^ t> b, enhancing the corresponding cross sections 
resulting to cosmologically acceptable relic densities as already discussed Such points 
allow for Ml/2, "t.o as large as ~ 450 GeV and stay comfortably well as far as the process 
6 — > s 7 is concerned, which is not in conflict with large tan (3 values as long as // > 
[28,43]. Since large values of tan/3 are compatible with Yukawa coupling unification, the 
previous discussion shows that the possibility of obtaining acceptable x relic densities in 
the coannihilation free region is feasible in such schemes. If Yukawa coupling unification 
is enforced, the input 6-quark pole mass should be lowered to vahics that are marginally 
consistent with the experimental data. This has as an effect the increase of the value of 
tan f3. In fact by lowering the input value M^°'^ = 5 GeV, we were able to get relic 
densities within the cosmologically allowed domain 0.08 <Q^hl < 0.22 for tan/3 ~ 50, 
without the need of invoking the coannihilation mechanism as is done in Ref. [39] . Note 
the important role the pseudoscalar Higgs boson plays in this case since it dominates 
the X X — > rf, bb reactions when the LSP's composition involves even a small Higgsino 
component. In figure 7, and in order to exhibit the tan /3 behaviour of the relic density. 



^■"^This requires the Axx coupling to be non-vanishing. This holds in regions of the parameter 
space where the LSP state has a non-vanishing Higgsino component. 
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we display a scattered plot of random points, for fixed Aq, M1/2 and random values of 
mo > 150 GeV, as function of tan P, for both signs of fi. All points displayed refer to the 
coannihilation free region under discussion. Actually for the /i > case, only a few points 
of the given sample are in the coannihilation region. We see that only a small number 
of points with < can marginally satisfy the cosmological constraints. However for 
> many such points exist for values of tan (3 which are around ~ 30. We recall that 
the bottom quark pole mass has been taken equal to 5 GeV which hardly allows for large 
values of tan /3. For this reason, and for the given sample, points beyond tan /3 ~ 35 (38) 
for // > (< 0) are absent in these figures. In the bottom figure, corresponding to // < 
case, we do not display points in the gap around tan (3 — 5, since we are close to a two 
hght Higgs threshold (see discussion in section III). 

EW precision data are in perfect agreement with the SM and hence also with su- 
persymmetric extensions of the SM which are characterized by a large supersymmetry 
breaking scale. In unconstrained SUSY scenarios the bounds put on sparticle masses 
from the EW precision data are not far from their lower experimental limits. In con- 
strained versions, such as the CMSSM which we study here, lower bounds on M1/2 can be 
established. In fact phenomenological studies of the weak mixing angle sin^^g// restrict 
Ml/2 to lie in the region M1/2 > 300 GeV if the combined small SLD and LEP data are 
used for sin^^g//^^- If in addition unification of gauge couplings at Mqvt is assumed then 
the lower bound is shifted to higher values (see Dedes et. al. in Ref. [14]), in the absence 
of high energy thresholds. Therefore in the context of the CMSSM it seems that EW 
precision data favour rather large M1/2 values in which case we are closer to the SM limit 
of Supersymmetry. The higher the M1/2 value the lower the is, and better agreement 
with the experimental data is obtained. Adopting a lower bound of about 300 GeV, sug- 
gested by the above reasoning, can have a dramatic effect for the allowed domain which 
hes entirely in the coannihilation free region. For low tan /3 (< 10) the cosmologically 
allowed region is severely constrained almost predicting the values of the soft masses. In 
fact Ml/2 is forced to move within the rather tight hmits M1/2 ~ 300 GeV — 340 GeV, 



The SLD data alone leave more freedom by allowing for lower M1/2 values. On the contrary 
small LEP data favour large M1/2 values. 
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while at the same time mo ~ 80 GeV — 110 GeV. For higher values of tan/? (~ 20) the 
upper bound on mo is sifted upwards by about 20 GeV (see for instance figure 1). This 
situation is depicted in figure 8a where in the (M1/2, ""^o) plane the dark-shaded area 
marks the cosmologically allowed region for values Aq = GeV, tan/3 = 5. The coan- 
nihilation free region under discussion lies above the line labelled by ruf = 1.25 rrij^. In 
this figure it is seen that by enforcing a more relaxed lower bound, M1/2 > 200 GeV, not 
excluded by SLD data, a relatively large portion in the (M1/2, ttiq) plane is allowed which 
also overlaps with regions in which neither color nor charge are violated (marked as "No 
CCB"^^) [44]. However for M1/2 > 300 GeV the allowed region, in the coannihilation 
free domain, is shrunk to a small triangle. 

The previously discussed bounds on M1/2, mo affect the mass spectrum of supersym- 
metric particles. For M1/2 = 300 GcV - 340 GcV, mo = 80 - 110 GeV, \Ao\ < 1 TeV 
and values of tan/3 < 10, wc have found the following bounds on the masses of the LSP 
and the lighter of charginos, staus, stops and Higgs scalars: 



^LSP 


: 115 (116) GeV 


- 130 (133) GeV, 


mg 


: 210 (218) GeV 


- 241 (250) GeV , 


mf^ 


: 122 (130) GeV 


- 157 (158) GeV, 




: 401 (403) GeV 


- 667 (687) GeV , 


mho 


: 96 (87) GeV - 


125 (122) GeV. 



These refer to the case fj, > (/i < 0). 

In order to see how the bound put on SUSY breaking parameters, and hence the 
sparticle masses are affected, if the more stringent cosmological limits quoted in Ref. [8] 
are employed, in the figure 8b we have drawn the same situation as in figure 8a with 
Qy, /iq = 0.12 ± 0.04. One notices that the decrease of the upper bound on Qy. from 
0.22 to 0.16 washes out the allowed points in the coannihilation free region, if the lower 
bound Ml/2 > 300 GeV is enforced. 



'The alert reader may notice that the overlap between the "No CCB" allowed region and the 
coannihilation region is of measure zero, at least for M1/2 less than 500 GeV. This trend may 
be very suggestive in looking for the physically sound region in parameter space. 
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Within the CMSSM the only option to evade the stringent bounds put on supersym- 
metry breaking parameters, and hence on sparticle masses, remains either to move to the 
large tan f3 regime, which we discussed previously, or to go to the coannihilation region 

< rrif^ < 1.25 in which case M1/2 is not actually bounded [31]. 

In the coannihilation region (ii) our results concerning the neutralino relic density 
do not hold any more. However the conclusions of our analysis and that presented in 
Ref. [31] can be both combined to infer information on the actual relic density, ftx^o^ 
from the one have calculated which we shall hereafter denote by Q| h^. Using the findings 
of this reference we can express the actual relic density as 

= R{AM) Q| , (41) 

where the reduction factor R{AM) depends on AM = (m^^ — m^)/m^ and is plotted in 
figure 9. It is seen that R{AM) smoothly interpolates between ^ 0.1 and 1.0 for values 
of AM in the range 0.00 — 0.25. The above equation is a handy device and reproduces 
the results cited in Ref. [31]. The cosmologically allowed domain shown in figures 8 
has been actually drawn using this equation. In figure 10 we see how the contours of 
figure 1 are distorted when Eq. (41) is implemented. Notice the change of the shape at 
the bottom of the figure where the mass of tr starts approaching that of the LSP. 

In the scattered plot of figure 11 we show all points of the random sample which were 
previously used for the production of figure 7, which lie strictly within the coannihilation 
region. These points were not displayed in the figure 7. In the figures at the top the 
vertical axis refers to values of the relic density which is based on our own calculation 
(0- hi). The second set of the figures, at the bottom, shows how some of these points 
collapse, if the Eq. (41) is used, falling within the cosmologically allowed stripe allowing 
for high Ml/2 (and mo), values. The vertical axis now refers to the actual relic density 

V. CONCLUSIONS 

In this paper we have evaluated the relic neutralino abundance in view of recent 
cosmological data which support evidence for a flat and accelerating Universe. The 
acceleration is mainly driven by a non-vanishing cosmological constant which weighs 
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about 2/3 of the total matter-energy density of the Universe. Such a large contribution of 
the cosmological constant (vacuum energy) pushes the matter density, and consequently 
the CDM density, to relatively small values f2cDM — 0-15 ± 0.07, constraining the 
theoretical predictions of supersymmetric extensions of the SM model. 

Supersymmetric theories, with i?-parity conservation, offer a comprehensive theoret- 
ical framework which provide us with a good candidate for the Dark Matter particle, the 
LSP, which turns out to be the hghtest of the neutralinos. The bound Q^^^ hi ~ 0.15±0.07 
shows preference towards low values for the effective supersymmetry breaking scale 
-^susY, which in conjunction with electroweak precision measurements, pointing to the 
opposite direction favoring rather large values for Msusy, put severe constraints affecting 
supersymmetric predictions. 

We have undertaken the calculation of the relic density in the context of the CMSSM, 
with radiatively induced breaking of the electroweak symmetry and universal boundary 
conditions for the soft supersymmetry breaking parameters in which the LSP plays the 
role of the Dark Matter particle. 

Our analysis have revealed the following: 

Although the cosmological data do not rule out corridors in the (mo, M1/2) plane in 
which the LSP is light, with substantial Higgsino mixing, with no bound put on sfermion 
masses, nevertheless such regions be excluded in view of the latest experimental data 
from chargino searches. 

Towards the large M1/2 regime we have found that in the cosmologically interesting 
domain, M1/2 cannot exceed ^ 340 GeV, while at the same time mo < 200 GeV. 
These bounds are obtained provided one stays within the region 1.25 rriy. < nif^ where 
coannihilation processes do not play any significant role. Putting a lower bound on M1/2 
suggested by EW precision data can have a dramatic effect on the allowed mo, M1/2 
values. If for instance, based on phenomenological studies of the electroweak mixing 
angle, we impose M1/2 > 300 GeV then mo, M1/2 are restricted to he within the tight 
hmits Ml/2 ^ 300 GeV -340 GeV, mo ~ 80 GeV -130 GeV. These limits are insensitive 
to the choice of the parameter Aq and hold as long as tan /3 < 30. If, as other analyses 
suggest, the more restrictive cosmological data are imposed, il^ = 0.12 ± 0.04, then 
there are no allowed points in the region m^^ > 1.25 m^j. for M1/2 > 300 GeV. 
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Within CMSSM there are two ways to reconcile the experimental information from 
EW and cosmological data with values of mo and M1/2 that he outside the strict bounds 
quoted above. We have either to go to the large tan f3 (with n > 0) regime, while staying 
within 1.25 < rrif^, or move inside the narrow band < < 1.25m^ in which 
case tr, the next to LSP sparticle, is almost degenerate in mass with the LSP and tr — x 
coannihilation processes are relevant to keep neutralinos in equilibrium. 

In the first case the pseudoscalar Higgs boson A plays an essential role. Depending on 
the inputs its mass may approach 2111^ while the Arf, Abb couplings are large as being 
proportional to tan (3. Both effects make the pseudoscalar Higgs exchange dominate the 
reactions xx tt", bb, and enhance the corresponding cross sections, resulting to relic 
densities which are compatible with the cosmological data. It is worth pointing out that 
large tan/? values, for fi > 0, are compatible with the CLEO data for the process b 57. 
In addition since large values of tan f3 are compatible with Yukawa coupling unification, 
this mechanism may offer the possibility of obtaining cosmologically acceptable x "celic 
densities in the coannihilation free region 1.25 < mf^, in such unification schemes. 

The second possibility to make the recent astrophysical data compatible with values 
of mo, Ml/2 outside the narrow domain quoted above, is to move to the coannihilation 
region < mf^ < 1.25 m^. In this region the tr — x coannihilation effects enhance 
(avrei), lowering significantly the values of the neutralino relic density. By using the 
findings of Ref. [31] we have found a handy way to relate the actual relic density fly. to 
that calculated using the traditional way, Q^, in which coannihilation reactions are not 
counted for. We find that in the region of the parameter space in which the LSP is nearly 
degenerate with the next to the LSP particle, namely Tr, no upper limit is imposed on 
the parameter M1/2. Given a value for M1/2 the parameter mo is however constrained 
to lie within a narrow band which is dictated by < mf^ < 1.25 m^^. 
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TABLES 



TABLE L Sample of values for the ratio gg/^ Qo for an LSP mass equal to 100 GeV. The 
masses of the remaining particles are as described in the main text. 



X 


go/2 ql 


.08 


2.47 X 10^ 


.07 


1.62 X 10^ 


.06 


1.95 X 10^0 


.05 


0.62 X 10^2 


.04 


1.05 X 10^^ 


.03 


0.51 X 10^8 
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Appendix: Supersymmetric conventions 



The supersymmetric Lagrangian we are using in this paper has a superpotential given 

by 

W ^ htH^ eQU" + h e Q D" + hr e L E" + t Hi (A.l) 

where the elements of the antisymmetric 2x2 matrix e are given by ei2 = — ei2 = 1 . 
In the superpotential above we have only shown the dominant Yukawa terms of the third 
generation. 

The scalar soft part of the Lagrangian is given by 

^scalar — ^ ] '^i 

i 

- {At ht eQU" + Ai,hb e Q D'= + K e L E" + h.c) 
+ (mlHj eHi + h.c) , (A.2) 

where the index i in the sum in the equation above runs over all scalar fields and all 
fields appearing denote scalar parts of the supermultiplets involved. 
The gaugino fields soft mass terms are given by 

jC-gaugino ^ { B B + + M3 (5 G + h.c). (A.3) 

In this equation B, W^'^\ G are the gauge fermions corresponding to the SU (2), 

and SU (3) gauge groups. 

For comparison with other notations [27,45] it is perhaps useful to remark that co- 
variant derivatives in this paper are defined by 

D,^d,-ig T(^) Af . 

Thus there is a sign difference in the gauge couplings used in this paper and in 
Refs. [27,45]. Besides that, the gaugino fields we use through differ in sign from those 
used in those papers and the parameters Mj and n are opposite in sign too. These 
remarks set the rules of passing from one notation to the other. 
In the B, W^^\ iH^, iH2, basis the neutrahno mass matrix is 
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M 



N 



(A.4) 



/Ml g'vl/^/2 -g'v2/V2\ 

M2 -gvi/V2 gv2/V2 

g'vi/V2 -gvi/V2 -ji 

\-g'v2/V2 gv2/V2 -/^ 

In this expression the tangent of the angle (3 sets the ratio of the v.e.v's of the two 
Higgses tan/3 = {H^) / (H,) . 

The mass eigenstates (x? 2,3,4) neutrahno mass matrix A4 n are written as 



/ 



O 







1 ^ \ 














\xlj 




\ iH'2 I 



(A.5) 



and 



M. nO — Diag {rn^o , m^o , m^o , m^oj , 



(A.6) 



where (9 is a real orthogonal matrix. Note that when electroweak breaking effects are 
ignored O can get the form 



o = 



0, 



02 j_ f 

V2 / 



(A.7) 



The chargino mass matrix can be obtained from the following Lagrangian mass terms 

(A.8) 



rmass 
charginos 



where we have defined = ^^^^Ti^^^^ and 



M 



c 



M2 - g V2 
- gvi 1^ 

Diagonalization of this matrix gives 

m^, 



(A.9) 



UMcV^ = 



"Xl 

m 



X2 



(A.IO) 



Thus, 



^charginos ~ I^XiXlXl ^X2X2X2 



(A.ll) 
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The Dirac chargino states Xi,2 are given by 

The two component Weyl spinors are related to W^, iH^, iH^ by 

K£)KaO ■ (*'-.'*r)c*MAr.Aj). (a.i3) 

The gauge interactions of charginos and neutrahnos can be read from the following 
Lagrangian^^ 



- 9 + W-J^^) + eA,J^^^ + -Z^Jf . (A.14) 

sc 



In the equation above s = sin6'vF, c = cos6'iy . Also, 



(A.15) 



The currents J^, and are given by 

^ £07'^ [VlV^, + VnVSl X^ a = 1...4, ^ = 1, 2 , (A.16) 
where Vl,r — and 

V^, = ~^0saU:,-02aU:,. (A.17) 

The electromagnetic current J^^ is 

JL = li7'^Xi + ^27^^X2 . (A.18) 
Finally, the neutral current can be read from 

J'z ^ ^iY [VlA^ + VnJ^^ Xi + [PLB^ab + -PrB^i] , (A.19) 
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In our notation e =eIectron's charge. 



33 



with 

= -B'a, . (A.20) 

Note that since = ~^ab neutrahno contribution to can be cast into the form 

J'z = -^^a^KTVxg) ■ (A.21) 

For the calculation of the X X ~^ f f cross sections we need know the chargino and 
neutralino couplings to fermions and sfermions. The relevant chargino couplings are 
given by the following Lagrangian terms 

>C = ^ ^- {Vl 4^ + Vr b0 ff;+t h in a{f + Vr b{f) f ~f* + {h.c) . (A.22) 

In this, Xi = 1, 2) are the positively charged charginos and Xi the corresponding charge 
conjugate states having opposite charge. / , /' are "up" and "down" fermions, quarks 
or leptons, while fi , f[ are the corresponding sfermion mass eigenstates. The left and 
right-handed couplings appearing above are given by 

In the equation above hf , hf/ are the Yukawa couplings of the up and down fermions 
respectively. The matrices K-^'-^' which diagonalize the sfermion mass matrices become 
the unit matrices in the absence of left-right sfermion mixings. For the electron and 
muon family the lepton masses arc taken to be vanishing in the case that mixings do not 
occur. In addition the right-handed couplings, are zero. 
The corresponding neutralino couplings are given by 

>C = ^ X° {Vl ai/ + VRbi^)ff;+t x° {Vl af/' + Vr 6^/) /' f]* + {h.c) . (A.23) 

The left and right-handed couplings for the up fermions, sfermions are given by 
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while those for the down fermions and sfermions are given by 
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FIG. 1. The LSP relic density in the (mo,Mi/2) pi ane for given values of 

tan/3 and sign of /i. Grey tone regions, from darker to lighter, designate areas 
in which the LSP relic density takes values in the intervals: 0.00 — 0.08, 0.08 — 0.22, 
0.22 - 0.35, 0.35 - 0.60 and 0.60 - 1.00 respectively. In the blanc area > 1.0. 

The area marked by "TH" is theoretically excluded (see main text). The area 
labelled by < 95 GeV is excluded by chargino searches. Crosses denote points 
for which thresholds or poles are encountered. 
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FIG. 2. The relic density as function of mo for fixed values of the remaining 
parameters. The solid, dashed and dot-dashed lines correspond to M1/2 = 170, 200 
and 400 GeV respectively. 
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FIG. 3. The relic density as function of ^o- The hnes are as in figure 2. 
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FIG. 4. The relic density as function of tan (3. The hnes are as in figure 2. 
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FIG. 5. The relic density as function of Mi/2- Crosses denote points that are near 
thresholds or poles. The solid (dashed) line corresponds to mo = 150 (200) GeV. 
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FIG. 6. Scattered plot of the relic density versus M1/2 from a sample of 4000 
random points in the parameter space. Low Mx/2 values are excluded by chargino 
searches. All points shown are in the coannihilation free region. Only the points with 
relic density less than 1.5 are shown. The grey tone region within the cosmologically 
allowed stripe designates the region which agrees with EW precision data (see main 
text). The horizontal dashed lines mark the limits 0.08 < J^j^^o ^ 0.22 . 
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FIG. 7. Scattered plots of the relic density versus tan (3 from a sample of random 
points with fixed = GeV, M1/2 = 300 GeV. The points shown fall within the 
coannihilation free region. The two horizontal dashed lines, as in figure 6, mark the 
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FIG. 9. The reduction factor R{AM) as function of AM = (mf^ — m^)/m^. 



48 



A =0, tanj8=5, 



A =0, tanj8=5, 



/2<0 




iliimliltliMNNlliltliNiiNfiiiilrl^Mi 



TH 



100 



200 
^1/2 



'I 

300 400 

(GeV) 



' I ' ' 
500 600 



600 



500 



.400 



> 
O 



300 



200 



100 




A =0, tanj8=20, /i,>0 



_L 



> 
un 



_L 



_L 



\ X X 

: X X 

XX 
., X X 
X X 
X X 

S,>< :x: : 
X x 
X X 

siSjiiSiiii 



X 
X 

m 

TH 



' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' 
100 200 300 400 500 600 

M,/2 (GeV) 



> 
O 



600 
500 
400 
300 
200 
100 




> 

a 
un 



tanj8=20, /i,<0 

_J I I \ I I I I \ I I I L 



X X 



XX 
X X 



X X 



X X 

X X 

X X 

X X 



X 

ills 



TH 



' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' 

100 200 300 400 500 600 

M,/2 (GeV) 



FIG. 10. The LSP relic density Q^h^ in the (mo,Mi/2) pl ane for eiven values 
of tan/3 and sign of fi when coannihilation effects are taken into account. The 
inputs are the same as in figure 1. 
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FIG. 11. The first two figures on the top show the points, from a random sample 
of 4000 points, which lie entirely within the coannihilation region. The vertical 
axis refers to 0| /^q (see main text). The figures at the bottom represent the same 
situation for the actual relic density h^. 
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